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Impact  signals  are  decomposed  on  orthogonal  wavelet  bases  or  wavelet  packet  bases. 
Based  on  these  orthogonal  decompositions,  the  energy  distribution  of  an  impact  signal 
can  be  defined  in  two  ways:  (a)  with  respect  to  the  natural  order  of  nodes  or  the 
frequency  index,  representing  the  energy  distribution  at  each  node  or  frequency  index;  (b) 
with  respect  to  the  time  position  and  the  frequency  index,  providing  an  energy  map  over  a 
time-frequency  plane.  By  comparing  the  energy  distributions  between  the  signals  from 
tests  and  simulations,  finite  element  impact  models  can  be  validated.  A  finite  element 
automobile  crashworthiness  model  is  used  to  illustrate  this  procedure  for  the  model 
validation. 


INTRODUCTION 

For  a  stationary  random  signal,  the  autospectral  density  function  (power  spectrum)  describes  its  energy  distribution 
with  respect  to  frequency.  For  a  non-stationary  random  signal,  the  energy  distribution  with  respect  to  frequency  is 
time-dependent.  Thus,  the  energy  distribution  is  a  function  of  time  and  frequency,  which  can  be  viewed  by  a  map 
over  the  time-frequency  plane.  There  are  several  methods  for  decomposing  non-stationary  signals  and  generating 
the  energy  distribution  map  [1].  Among  these  are:  (a)  the  short-time  Fourier  transform  (STFT),  which  generates  a 
spectrogram;  (b)  the  Wigner-Ville  method,  creating  a  time-frequency  distribution;  and  (c)  the  wavelet  analysis, 
providing  a  time-scale  distribution  [2]. 

Impact  signals,  such  as  automobile  crash  or  aircrew  ejection  responses,  have  several  unique  features.  They  are 
transient,  as  an  impact  event  usually  takes  place  in  short  time  duration.  The  dominant  or  major  component  of  an 
impact  response,  such  as  the  gross  motion  of  a  vehicle  in  a  frontal  impact,  can  be  considered  as  deterministic. 
Some  impulses  or  short-time  vibrations,  resulting  from  secondary  impact  or  structure  failure,  may  occur  in  certain 
time  spans.  The  impact  signals  are  usually  contaminated  with  random  noise.  Whereas  an  impact  signal  can  be 
broadly  considered  as  a  non-stationary  signal,  a  new  method  for  energy  distribution  analysis  is  desired  in  order  to 
describe  its  unique  features. 

As  a  new  tool  for  signal  analysis,  wavelets  are  localized  in  both  frequency  and  time  domains,  which  match  the 
major  characteristics  of  impact  responses.  Therefore  wavelets  are  used  in  this  paper.  Impact  responses  will  be 
decomposed  on  orthogonal  wavelet  bases  or  wavelet  packet  bases.  Based  on  these  orthogonal  decompositions,  the 
energy  distribution  of  an  impact  signal  will  be  defined  and  the  method  for  the  energy  distribution  analysis  will  be 
developed.  The  energy  distribution  analysis  will  be  performed  on  automobile  crash  signals  from  an  actual  test  and 
from  computer  simulations.  As  a  practical  application,  the  energy  distribution  will  be  used  in  the  validation  of 
finite  element  automobile  crashworthiness  modeling. 

WAVELET  AND  WAVELET  PACKET  DECOMPOSITION 
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In  a  wavelet  basis,  a  signal  S(t)  can  be  decomposed  as  [3-5] 


iS  —  +  Z)j 

=  A2  H"  D2  "b 

j^J 


(1) 


where 

(2) 

k 

and 

Dj{t)=Y.^jk¥jit-k).  (3) 

k 

Here,  Aj  and  D j  are  referred  to  as  the  approximation  and  the  detail  of  S  at  level  j ;  and  y/ j{t)  are  the 

scaling  function  and  the  wavelet  function  at  level  j  for  reconstruction;  and  Cjj^  and  d ,  given  by  wavelet 

transforms,  are  scaling  function  coefficients  and  wavelet  coefficients  at  level  j  and  time  shift  k ,  respectively.  The 
decomposition  of  a  signal  in  a  wavelet  basis  is  illustrated  in  Fig.  1(a). 


nn-i 


Ai 

A7- 

D2 

A3 

ih 

Di 


(a)  Wavelet 


(b)  Wavelet  packet 


Fig.  1 .  Decomposition  in  a  wavelet  basis  and  a  wavelet  packet  basis 


The  total  band  is  divided  unevenly  using  a  wavelet  analysis.  To  obtain  the  decompositions  over  evenly  spanned  sub¬ 
bands,  the  wavelet  packet  analysis  can  be  used  [3,4].  As  illustrated  in  Fig.  1(b),  the  wavelet  packet  analysis  is  a 
generalization  of  the  wavelet  analysis.  In  the  orthogonal  wavelet  decomposition  procedure,  successive  approximation 
coefficients  are  split  into  two  parts,  but  successive  details  are  never  re-analyzed.  In  the  corresponding  wavelet  packet 
situation,  each  detail  coefficient  vector  is  also  decomposed  into  two  parts  using  the  same  approach  as  in 
approximation  vector  splitting.  This  offers  a  more  complete  analysis  of  a  signal. 


In  a  wavelet  packet  basis,  the  decomposition  of  the  signal  S(t)  at  level  j  can  be  expressed  as 


2^-1 


n=0  k 


(4) 


2 


where  VVy^(0  are  wavelet  packet  atoms  [3,4];  are  wavelet  packet  coefficients;  n  is  the  frequency  index  which 

is  in  accordance  with  the  natural  order  of  the  nodes  at  this  level;  and  k ,  the  time  shift,  is  the  position  index.  For  a 

2^-1 

J  -level  decomposition,  from  a  complete  wavelet  packet  decomposition  tree  there  are  more  than  2  different 
ways  to  encode  a  signal  [6],  for  which  a  more  general  expression  is 

^  ^  ^  9  jnk  ^jn  (^  “  ^)  »  (5) 

U.n)el  k 

where  I  is  an  index  set  which  is  chosen  as 


(6) 


such  that  intervals  [2*^'  ,2-^'  + 1))  are  disjointed  and  cover  the  entire  interval  [0,oo)  [7]: 


U[2^'/2,.,2^'(«,.+1))  =  [0,co). 

/=0 


(7) 


One  such  decomposition,  which  is  picked  from  the  decomposition  tree  of  Fig.  1  (b),  can  be,  for  instance, 

S  =  ‘^‘^2  ^^^3  ^^^3  ^^2  •  (^) 

As  the  number  of  ways  for  the  decomposition  of  a  signal  in  wavelet  packet  bases  may  be  very  large  and  since  explicit 
enumeration  is  generally  unmanageable,  it  is  necessary  to  find  the  optimal  decomposition  in  terms  of  certain  entropy 
criteria,  computable  by  an  efficient  algorithm  [3,4].  These  criteria  include  Shannon,  Threshold,  Norm,  and  Log 
Energy,  which  are  available  in  the  MATLAB  wavelet  toolbox  [4]. 

SIGNAL  ENERGY  DISTRIBUTIONS 

In  a  wavelet  analysis,  if  orthogonal  wavelets  are  used,  all  wavelets  should  be  orthogonal  to  the  scaling 

functions  ^(/  -  k)  .  Furthermore,  the  wavelets  y/ j  (/  -  k)  should  be  mutually  orthogonal,  and  the  scaling  functions 
^(/  -  k)  should  be  mutually  orthogonal  also  [3,  5].  This  leads  to  the  following  relations  [3]: 

(l)j{t-m)(l>j{t-n)dt  =  5{tn-n), 

(pj  it-m)i//j(t-n)dt  =  S(m-n),  (9) 

(t  - k)ii/j (J  - K)dt  =  5{j  -  J)d{k  - K) , 

where  S  is  the  Kronecker  delta  function.  Accordingly,  the  hierarchical  decomposition  of  Eq.  (1)  is  orthogonal.  That 
is, 

•  Aj  is  orthogonal  to  Dj ,  _2 , .  •  •  > 

•  Dy  is  orthogonal  to  for  j  ^  k  . 

It  follows  that  in  terms  of  the  energy  that  S(t)  contains. 
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=  U,  +  A 


=kir +iiAir+iiA 


(10) 


=Kir+ziia, 

j^J 

Furthermore,  from  the  relations  expressed  by  Eq.  (9), 

iisf=Z4+i;2;4. 

k  7=1  k 

which  describes  the  energy  distribution  of  a  signal  over  scales  (frequency  bands)  and  time  spans. 


(11) 


Wavelet  packet  functions  (atoms)  are  created  by  taking  linear  combinations  of  the  usual  wavelet  functions  [3,4,7]. 
The  wavelet  packet  bases  are  a  generalization  of  wavelet  bases  and  inherit  properties  such  as  orthogonality  and 
smoothness  from  their  corresponding  wavelet  functions  [7].  Therefore,  the  relations  similar  to  those  described  for  the 
wavelet  analysis  may  exist  in  the  wavelet  packet  analysis.  Among  these  is,  for  orthonormal  wavelet  packet  atoms. 


£  Wj„  it  -  it  -  l)dt  =  Sij  -  J)Sin  -  N)Sik  - 1) . 

If  the  decomposition  of  a  signal  is  performed  on  orthogonal  wavelet  packet  bases,  in  terms  of  the  signal  energy, 

„.„2 


=  ’ 


(12) 


(13) 


n=0  k 


which  corresponds  to  Eq.  (4),  and 


which  follows  from  Eq.  (5). 


(y,«)e/  k 


(14) 


The  signal  energy  distribution  with  respect  to  frequency  and  time  can  be  determined  using  Eq.  (13)  if  the  signal  is 
decomposed  at  a  certain  level  and  on  evenly  spanned  frequency  bands.  Likewise,  the  signal  energy  distribution  with 
respect  to  scale,  frequency,  and  time  can  be  analyzed  based  on  Eq.  (14)  when  the  signal  is  decomposed  on  an  optimal 
tree  that  may  spread  over  multiple  levels  and  on  unevenly  spanned  frequency  bands. 


ENERGY  DISTRIBUTION  ANALYSIS  FOR  MODEL  VALIDATON 


For  the  finite  element  (FE)  modeling  of  automobile  crashworthiness,  the  model  validation  is  often  conducted  in  terms 
of  entire  vehicle  motions  and  acceleration  responses  at  certain  critical  points  [8,9].  When  acceleration  responses  are 
compared  in  the  validation,  the  time  histories  measured  in  actual  tests  and  obtained  from  simulations  are  often  plotted 
curve  by  curve.  The  agreement  or  discrepancy  between  the  test  data  and  simulation  results  can  be  evaluated  by  visual 
inspection  of  the  amplitudes,  phases  or  timing  of  peaks,  and  pulse  shapes  of  both  signals.  While  visual  inspection  is 
intuitive,  direct,  and  easy,  it  can  be  qualitative,  subjective,  and  coarse.  Quantitative  and  objective  evaluations  based 
on  statistical  or  other  scientific  analyses  can  be  more  accurate  and  meaningful,  and  thus  are  desirable. 

The  energy  distribution  describes  the  variation  of  amplitude  with  respect  to  frequency  and  time.  The  comparison  of 
the  energy  distributions  of  two  signals  provides  a  description  of  the  agreement  and  discrepancy  between  them;  thus  it 
can  be  used  for  the  validation  of  a  computational  model  where  the  simulated  responses  using  the  model  are  compared 
with  actual  responses  from  tests.  When  a  pair  of  signals  from  test  and  simulation  are  decomposed  in  the  same  wavelet 
packet  basis  with  the  same  entropy  criterion,  the  best  decomposition  trees  for  them  may  be  different.  This  means  that 
the  signal  energy  distributions  based  on  respective  best  decomposition  trees  may  not  be  comparable.  Therefore,  both 
signals  need  to  be  decomposed  on  the  same  tree,  in  order  for  their  energy  distributions  to  be  comparable.  As  a  choice, 
both  signals  can  be  decomposed  at  the  same  level  in  a  wavelet  basis  (Eq.  (1))  or  a  wavelet  packet  basis  (Eq.  (4)). 
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Fig.  6.  Energy  distributions  with  respect  to  k 
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(d)  Right-rear  cross  member 
Fig.  4.  Energy  distribution  with  respect  to  n  and  k 
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Fig.  5.  Energy  distributions  with  respect  to  n 
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Energy  Distribution  Analysis  and  Model  Validation 

The  energy  distributions  of  these  signals  with  respect  to  the  natural  order  and  time  location,  which  are  calculated 
using  Eq.  (16),  are  displayed  in  Fig.  4.  This  figure  provides  a  three-dimensional  view  of  the  signal  energy 
distributions.  For  the  acceleration  responses  at  the  locations  of  engine  bottom  (Fig.  4(a))  and  engine  top  (Fig.  4(b)), 
the  energies  of  both  signals  from  test  and  simulation  concentrate  on  the  first  node  with  natural  order  a?  =  0 
corresponding  to  the  lowest  frequency  band,  and  spread  over  th*e  time  span  from  the  position  k  =  5  to  k  =  \5  .  For 
the  acceleration  responses  at  the  locations  of  the  left-  and  right-rear  cross  member  (Fig.  4(c)  and  (d)),  both  signal 
energies  are  primarily  contained  at  the  first  node,  but  there  are  significant  amounts  of  energy  distributed  over  high 
order  nodes,  especially  for  the  response  at  the  left-rear  cross  member  from  the  simulation.  The  energies  at  these 
nodes  spread  over  a  wide  time  span. 

To  find  the  energy  distributions  at  each  node,  the  energy  distributions  with  respect  to  the  natural  order  are  calculated 
using  Eq.  (15)  and  are  illustrated  in  Fig.  5.  Since  the  energy  distributions  over  high  order  nodes  in>32)  are 
insignificant,  as  shown  in  Fig.  4,  the  energy  distributions  in  Fig.  5  are  displayed  only  for  the  natural  order  up  to 
A?  =  3 1 .  The  entire  range  is  from  ai  =  0  to  a?  =  63  .  In  addition  to  providing  a  more  clear  view  of  the  energy 
distributions  at  each  node  as  compared  to  Fig.  4,  Fig.  5  offers  the  comparison  between  the  signals  from  test  and 
simulation.  As  far  as  the  principal  component  of  energy  that  is  placed  on  the  first  node  is  concerned,  the  agreement 
between  the  test  data  and  simulation  results  using  the  FE  model  is  sound  for  the  acceleration  response  at  the  right- 
rear  cross  member,  good  at  the  engine  bottom,  fairly  good  at  the  engine  top,  and  still  close  at  the  left-rear  cross 
member  where  considerable  amounts  of  energy  are  contained  at  high  order  nodes  (Fig.  5  (c)). 

The  energy  contained  in  each  node  can  be  further  analyzed  from  the  perspective  of  its  distributions  over  time,  as 
shown  in  Fig  4.  To  obtain  a  more  intuitive  and  clear  view,  based  on  Eq.  (16),  two-dimensional  graphs  are  derived 
from  Fig.  4  for  the  nodes  that  have  major  signal  energy  distributions.  These  graphs  are  shown  in  Fig.  6.  As  far  as  the 
energy  contained  in  the  first  node  is  concerned,  the  simulation  with  the  FE  model  has  attained  a  sound  agreement 
with  Ae  actual  test  in  terms  of  the  amount  of  energy  in  this  frequency  band,  but  the  distributions  over  time  spans  are 
not  quite  agreeable.  For  other  higher-order  nodes,  these  distributions  are  basically  not  agreeable  between  the 
simulation  and  the  test. 
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CONCLUDING  REMARKS 


Based  on  the  decomposition  of  a  signal  on  wavelet  packet  bases,  two  types  of  signal  energy  distributions  have  been 
established:  (a)  the  energy  distributions  with  respect  to  the  natural  order  of  nodes  or  the  frequency  index  n ;  (b)  the 
energy  distributions  with  respect  to  the  time  position  k  and  the  frequency  index  n .  The  unique  characteristics  of  a 
wavelet  analysis  render  the  energy  distributions  based  on  wavelet  or  wavelet  packet  decompositions  to  be  especially 
suitable  for  the  analysis  of  impact  signals,  such  as  automobile  crash  responses.  The  gross  motion  component  of  an 
impact  signal  is  readily  represented  by  the  approximation.  The  short-time  impulses  and  vibrations  localized  in  both 
time  and  scale  (frequency)  can  be  expressed  easily  by  the  details.  If  necessary,  de-noising  and  compression,  two 
powerful  operations  provided  by  a  wavelet  analysis,  can  be  employed  to  eliminate  noise  and  insignificant 
components,  so  that  the  energy  distribution  can  focus  on  significant  components. 

It  is  important  to  realize  that  the  energy  distributions  over  the  time-frequency  plane  based  on  wavelet  or  wavelet 
packet  decompositions  may  look  quite  different  when  different  wavelets  are  used.  Unfortunately,  this  problem  occurs 
in  other  methods  (such  as  STFT)  also  [2].  Therefore,  caution  needs  to  be  taken  when  a  wavelet  time-frequency 
spectrogram  is  interpreted.  The  characteristics  of  the  wavelets  used  in  the  decomposition  need  to  be  taken  into 
consideration  also. 

Due  to  transient  and  non-stationary  characteristics  of  automobile  crash  responses,  validation  methods  and  criteria 
based  on  conventional  statistic  analysis,  correlation  analysis,  or  spectral  analysis  used  for  stationary  signals  may  be 
neither  appropriate  nor  efficient.  Comparing  energy  distributions  between  a  pair  of  signals  from  tests  and  simulations 
provides  a  way  to  evaluate  simulation  results  and  to  validate  the  model  from  certain  perspectives.  The  distributions  of 
type-b  provide  more  detail  and  complete  information  of  a  signal  than  the  distributions  of  type-a.  Accordingly,  the 
validation  based  on  the  energy  distributions  of  type-b  is  more  comprehensive  and  rigorous  than  the  validation  based 
on  the  distributions  of  type-a. 

Whereas  the  energy  distributions  defined  above  provide  a  common  basis  for  the  comparison  of  two  signals,  they  are 
mainly  concerned  with  the  amplitudes  of  signals.  It  is  sometimes  necessary  to  compare  the  shapes  of  pulses  and  the 
timing  of  peaks  between  two  responses.  This  problem  can  be  treated  from  the  perspective  of  the  correlative 
relationship  between  two  signals  [14,15]. 

The  validation  of  a  1997  Honda  Accord  FE  model  for  full  frontal  impact  has  been  performed  in  this  paper.  The 
impact  responses  from  a  simulation  and  an  actual  test  have  been  evaluated  and  compared.  The  results  show  that,  from 
the  perspectives  of  the  energy  distributions  over  frequency  bands,  the  agreement  between  the  simulation  results  using 
this  FE  model  and  the  test  data  is  good  for  the  gross  motions  at  the  locations  of  the  engine  bottom,  engine  top,  and 
right-rear  cross  member,  and  fairly  good  for  the  gross  motion  at  the  location  of  the  left-rear  cross  member. 
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